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We propose a nonadiabatic time-dependent spin-density functional theory (TDSDFT) approach 
for studying the single-electron excited states and the ultrafast response of systems with strong 
electron correlations. The correlations are described by the correlation part of the nonadiabatic 
exchange-correlation (XC) kernel, which is constructed by using some exact results for the Hubbard 
model of strongly correlated electrons. We demonstrate that the corresponding nonadiabatic XC 
kernel reproduces main features of the spectrum of the Hubbard dimer and infinite-dimensional 
Hubbard model, some of which are impossible to obtain within the adiabatic approach. The theory 
may be applied for DFT study of strongly correlated electron systems in- and out-of-equilibrium, in- 
cluding the important case of nanostructures, for which it leads to a dramatic reduction of necessary 
computational power. 
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/77.trorfMction. -Description of the electronic properties 
of the systems that contain localized d- and f-orbitals re- 
mains one of the most important problems in condensed- 
matter physics for both extended and finite systems. The 
extended systems form a large class of materials with 
many exotic properties with in turn potentially many 
technological applications. Examples include cuprate 
high-temperature superconductors, heavy fermion mate- 
rials and manganites. One can also expect many unusual 
properties in the much less explored cases of molecules 
and nanostructures, in which correlation effects may be 
even more enhanced owing to space constraint. More- 
over, in modern nanotechnological applications, in which 
the distance between the atoms can be tuned (for ex- 
ample, by putting the atoms on a substrate), even s- 
and p-electron systems can find themselves in a strongly- 
correlated regime when the interatomic separation is 
large enough. Such systems show many surprises, like 
recently observed unusual antiferromagnetic ordering in 
small Fe chains,- metal-insulator transition in Au and Fe 
chains'^ etc. The excited states and the nonequilibrium 
properties of such systems, including the ultrafast (femto- 
and atto-second) response, are extremely important from 
the point of view of the modern trend of "smaller and 
faster." Correct understanding of the nanosystems and 
molecules with strong electron correlations may also shed 
light on the general properties of strongly correlated ma- 
terials, including the local correlations and nonhomoge- 
neous order in extended systems. 

Most of the progress in studies of correlated systems 
has been made by using many-body approaches. Two of 
the most powerful of them are the Bethe ansatz in the 
ID case and the dynamical mean-field theory (DMFT) 
approach, the latter of which is exact in the limit of 
infinite dimensions'^ and which appears to be also a 
good approximation in the 2D and 3D cases (reviews of 
both approaches are given in Refs. Eland EJ correspond- 
ingly). These approaches have also been generalized to 



the nonequilibrium case,'^ allowing one to study the ex- 
citations and the nonlinear response of the correspond- 
ing systems. Since standard DFT approximations fail 
to describe properly most of the strongly-correlated ef- 
fects, nowadays DMFT combined with DFT is used for 
this purpose.EEl This method allows one to describe the 
spectral, optical and magnetic properties of bulk and 
layered materials (see, for example, reviews [lOllip . In 
DFT-hDMFT, aU properties of the "non-correlated" sys- 
tem (system geometry, bandstructure, etc.) are obtained 
within DFT (usually with the LDA or the GGA approx- 
imation), and the correlation effects are taken into ac- 
count by solving the corresponding effective Hubbard 
model. Recently, it has been shown that a similar ap- 
proach can be successfully applied to nanosystems .I^^Hlfil 
Unfortunately, the combined DFT-I-DMFT calculations 
are computationally demanding even for the equilibrium 
case, especially for finite systems, which have a number 
of non-equivalent atoms. 

For this reason it would be very useful to develop a 
TDDFT formalism with the corresponding XC potential 
that properly describes strongly correlated systems, in- 
cluding the excitations (in which case one needs to go 
beyond static DFT).-- Some progress in this direction 
has already been made. For ID systems, the correspond- 
ing adiabatic (static DFT) XC potentials were proposed 
in Refs. [^DH251 In particular, in the last two papers a 
DFT approach based on the Bethe ansatz has been pro- 
posed. Recently, a combined adiabatic TDDFT-DMFT 
approach applicable in the 3D case has been proposed, 
in Ref. In this approach, the XC potential is derived 
from the DMFT solution for the Hubbard model. The 
approaches were tested versus some exact solutions, and 
it was shown that they are successful in reproducing fea- 
tures of some of very important effects, including metal 
insulator transition and a temporal response of the sys- 
tems. Unfortunately, some of the effects turned out not 
to be properly described. In particular, as was shown in 
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Ref. the corresponding adiabatic potential fails to de- 
scribe correctly the response of a finite Hubbard system 
in strongly correlated regime when the system is close 
to half-filling and when the local Coulomb repulsion U 
is large. One might expect that the main shortcomings 
come from the adiabatic approximation. Indeed, as one 
notes, the adiabatic XC kernel fails even to reproduce 
the correct number of the excited states of the system. 
For example, in the case of the Hubbard dimer, it is 
known that the local Coulomb repulsion - in addition 
to the bonding-anti-bonding orbital excitation - leads to 
a satellite level at finite U (see, e.g., Ref. [25]) . This ex- 
tra state is related to the excited Hubbard "band" in 
the case of extended system (whence in fact comes the 
redistribution of the single-electron spectral weight to a 
higher energy (~ U). It is also known from the DMFT 
solution that the Hubbard systems demonstrate an extra 
quasi-particle spectral weight at zero energy (chemical 
potential).- Although all these states are important in a 
strongly correlated regime, none can be reproduced with 
the adiabatic TDDFT. Indeed, in the adiabatic case the 
solution of the Casida equatiorPS! leads to a shift of the 
single-electron levels , it does not contain new states. 

In this paper, we propose on the basis of some ex- 
act results for the Hubbard model a simple form of 
the nonadiabatic XC kernel which results in the single- 
electron spectrum of the Hubbard model that reproduces 
the main features of the spectra of both finite (dimer) and 
extended (infinite-dimensional) systems. This kernel can 
be easily implemented within the standard TDDFT codes 
for use in the strongly correlated case. 

The Hubbard dimer. -In order to obtain the XC ker- 
nel /_xCfT(T'(r, r', w) for the Hubbard dimer, we map the 
eigenvalue equation derived from the dimer Green's func- 
tion onto the corresponding TDDFT Casida eigenenergy 
equation.'2£l The Casida equation for the dimer has the 
following general form: 
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(the first part being the Hartree term, 
and /jfC(T,cr'(r,r',w) Fourier transform of 
SVxc<7[n]{'r,t)/Snc'{r' ,t') with respect to t-t'), 
eo = ~ £3 is the excitation energy of the free electron, 
eu,g and V'ii,g(r) are the corresponding bonding- and 
anti-bonding energies and wave functions. This equation 
has two solutions 
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where -I- corresponds to the singlet state, and — to the 
triplet one. Since the ground state of this system is sin- 
glet and the total spin of the isolated system is conserved, 
for definiteness we shall focus on the singlet state. 



The dimer Green's function can be found from the ex- 
act solution for the Hubbard model with the Hamilto- 
nian: 



(4) 



where c^^ and Cja are the creation and annihilation op- 
erator of electron with spin a on site i, rii^r is the corre- 
sponding number of particle operator, and t is the hop- 
ping parameter. One can find the exact single particle 
dimer's Green's function, which has the following form 
in the singlet (or, more generally, the "non-magnetic") 
case: 
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where the self-energies are 
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It is easy to show that, in this case, the corresponding 
eigenvalue problem has the following form: 



2 [/2 U^6t^^uy4_ 2 
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Hence, one can obtain an equation for the XC kernel for 
the Hubbard dimer: 

^nH + ^nM-^(i + ^;^), (9) 

where K^-a'i^) is defined in Eq. ([2]). Since we are inter- 
ested in the contribution to the energy due to correlation 
effects, we assume that the Hartree l/|r — r'|) and the 
exchange {K-^-]-) parts simply lead to the renormalization 
of the free-particle energy eg — > t in Eq. ([l]). 

In this case, one can arrive at the following separable 
form of the correlation portion of the local XC the kernel, 
which follows from Eq. (|9]): 



/ct,;(r,r',^)=5(r-r')^(^l + 
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were A = J / |V'gP(r)|i/;„p(r)drcir'. This kernel results 
in the exact excitation spectrum of the Hubbard dimei^^Sl 
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FIG. 1: Exact and DMFT spectra for the Hubbard dimer as 
a function of U/t. 



The spectra are presented in Fig. la. As follows from 
this Figure, there are two extra states Hubbard satellite 
peaks with energy ±3t at small U's. These states ap- 
pear only at finite U, and come from the redistribution 
of single-electron spectral weight. It is easy to see that 
one cannot obtain these new states from Eq. ^ with the 
static XC kernel. 

In the case of extended system, considered below, the 
DMFT approximation ~ 5ijY,{uj) is valid. Even 

though in the case of the dimer this approximation is 
not sufficiently accurate, we present the corresponding 
expression for the kernel for the discussion below: 
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As follows from the last two equations, they have the 
same expression at the high-frequency limit, but as w — >■ 
they differ. The DMFT solution has an extra energy 
peak around zero energy, which is characteristic for this 
approach (Fig.l). 

Extended systems.-ln these cases, the DMFT with the 
local self-energy approximation can be applied. We shall 



use the following exact result for the high-frequency elec- 
tron self-energy for the Hubbard model: 



(13) 



(see, for example Refs. I27l28jl . In the homogeneous case, 
one can construct the XC kernel by mapping the eigen- 
value equation 



or 
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onto the corresponding Casida equation. As in the dimer 
case one can find the equation which connects the XC 
kernel with the self-energy: 



K^^iuj) ^ i?e[E(w)]. 
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A more straightforward way to find the expression for 
the XC kernel is to compare the TDDFT and DMFT 
correlation energies: 
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Equations (13 1, (171 and ( |18[ ), together with the result 
for the dimer, Eq. ( |Io| ), can be used to construct the 
following " universal" (DMFT) function for the exchange 
part of the XC kernel for the extended systems: 



B2 



where i^[7io](i") is a functional of the ground state density 
no(r) (in general, of the spin parts) and i? ~ 3i in the 
case of dimer, while one can choose to be equal mean 
square of the kinetic energy , or more generally — 
asj., where a ~ 1. is average per site number of 
the spin-down electrons. It is important to note that 
we do not include the contribution of the static Hartree 
term Un^ to the self-energy, since it is canceled by the 
chemical potential at half-filling. With this XC kernel 
one can reproduce the main features of the spectrum of 
the infinite- dimensional Hubbard model: splitted (by 
energy U) bands and the zero energy quasi-particle peak, 
which disappears as U increases (Fig. 2). 

Conclusions.-As it is known from the DFT+DMFT 
study of real materials, the main effects of strong elec- 
tron correlations can be described by local, temporally 
resolved electron interaction. Therefore, one can expect 
that a spatially local, nonadiabatic XC kernel is suffi- 
cient for description of the strongly correlated systems 
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FIG. 2: a) The DOS for the infinite-dimensional 

Hubabrd model with the hypercubic DOS A{e) = 
{l/^/wt*)exp{—e^ /t*^) (t* = 2dt is the renormalized hop- 
ping, d - dimensionality of the system) and F[no]{r) — l, A=l, 
B — e^; b) the same in the low frequency case. 



in the TDSDFT case. From our analysis above it follows 
that one can expect the following nonadiabatic XC ker- 
nel can be used to describe the main effects due to strong 
electron-electron correlations within TDDFT: 



where B is a parameter proportional to the kinetic energy 
(B=3t in the case of dimer). The functional F[no-o](r) is 
defined by the static (" non-correlated" , DFT) spin den- 
sities. In order to model strong on-site correlations it can 
be chosen (for spin- up electron) in the form F[na-o](r) — 
n-fo(r)/(n-fo(r) -f- n^o(r)) (equal 1/2 in the case of the 
dimer). Indeed , this function has the maximum when 
the corresponding site (space point with a high charge 
density) is occupied by spin-up electron only. The fre- 
quency power a in the numerator is expected to be be- 
tween zero (in small systems, like dimer), and 1, as in the 
extended (DMFT) case. It is interesting to note that the 
case a=0 is in agreement with the exact asymptotic be- 
havior of the XC kernel /xc('^ — ^ oo) — >■ A+B/uj'^, where 
A and B are constants.- The most essential features of 
the proposed XC kernel are spatial locality, proportion- 
ality to the local Coulomb repulsion, and the oscillating 
in time interaction of the electrons with opposite spin, 
^ cxp{±iBt), where the frequency of oscillations B is 
proportional to the hopping (kinetic) energy. 



While further tests of the validity of the methodol- 
ogy presented here for larger scale systems, both in- 
and out-of-equilibrium is necessary, is on-going, the fre- 
quency dependence of the XC kernel in Eq. ( 20 ) is suffi- 



cient to describe the main features of single-orbital corre- 
lated systems. We expect the present approach to open 
the possibility of describing strongly correlated mate- 
rials within the standard TDDFT framework, extend- 
ing it beyond the metallic and semiconductor structures 
with plasmons, excitons and other excitations (see, e.g., 
Refs.lSHEIl)- 
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